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Abstract 

The complex structure of a typical stratus cloud base height (or profile) 
time series is analyzed with respect to the variability of its fluctuations and 
their correlations at all experimentally observed temporal scales. Due to the 
underlying processes that create these time series, they are expected to have 
multiscaling properties. For obtaining reliable measures of these scaling prop- 
erties, different methods of statistical analysis are used herein : power spectral 
density, detrended fluctuation analysis, and multifractal analysis. This broad 
set of diagnostic techniques is applied to a typical stratus cloud base height 
(CBH) data set; data were obtained from the Southern Great Plains site of 
the Atmospheric Radiation Measurement Program of the Department of En- 
ergy from a Belfort Laser Ceilometer. First, we demonstrate that this CBH 
time series is a nonstationary signal with stationary increments. Further, two 
scaling regimes are found, although the characteristic laws are quite similar 
ones. Next, the multi-affine scaling properties are confirmed. The scaling 
properties of the cloud base height profile of such a continental stratus are 
found to be similar to those of the marine cloud base height profiles studied by 
us previously. Some physical interpretation in terms of anomalous diffusion 
(or fractional random walk) is given for the continental case. 

Keywords : stratus cloud, cloud base height, fluctuations, correlations, power spectrum, 
detrended fluctuation analysis, multifractals 
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I. INTRODUCTION 



The state of the atmosphere is governed by the classical laws of fluid motion and exhibits 
a great deal of correlation at various spatial and temporal scales. Knowing these space- 
and time-scale-dependent correlations is crucial, for example, in order to understand the 
short and long term trends in climate. In particular, clouds play an important role in the 
atmospheric energy budget. In order to model and predict climate successfully, we must be 
able to both describe the effects of clouds in the current climate and predict the complex 
chain of events that might modify the distributions and properties of clouds in an altered 
climate. Moreover, in order to improve the parameterization of clouds in climate models, it 
is important to understand the cloud properties and their changes within the cloud. 

Modeling the impact of clouds is difficult because of their complex shapes, spatial dis- 
tributions and varying particle size distributions and because of their differing effects on 
weather and climate. Clouds can reflect incoming sunlight, and so contribute to cooling, 
but they also absorb infrared radiation leaving the earth, and so contribute to warming. 
High cirrus clouds, for example, may have a nonnegligible impact on atmospheric warming. 
Low-lying stratus clouds, which are frequently found over oceans, can contribute to warming 
as well.0 

A variety of physical processes takes place in the atmospheric boundary layer. At time 
scales of less than one day, significant fluxes of heat, water vapor and momentum occur 
due to entrainment, radiative transfer, and/or turbulence. 0i The turbulent character of the 
motion in the atmospheric boundary layer (ABL) is one of its most important features. 
The turbulence^ can be caused by a variety of processes, among them thermal convection 
and mechanical generation by wind shear.lfii This complexity of physical processes and 
interactions between them creates a variety of atmospheric responses. In particular, in a 
cloudy ABL, the radiative fluxes produce local sources of heating or cooling within the mixed 
layer and therefore can greatly influence its turbulent structure and dynamics. Moreover, 
variations in the turbulent structure and dynamics of the clouds cause subsequent changes 
in the cloud boundaries, especially in the height of cloud base. These variations lead to a 
complex structure for the cloud base height data (time) series. To analyze different aspects 
of the variability of its fluctuations and correlations at all experimentally observed temporal 
and spatial scales, one needs to apply a variety of diverse techniques of statistical analysis 
to the retrieved time series data. In what follows, we briefly present a few of these methods: 
a traditional technique - power spectral density, and rather new techniques, like a detrended 
fluctuation analysis, and multifractal analysis and then apply them to cloud base height 
data. 

II. SCALING 

A. Power spectral density 

The power spectral density S(f) of a signal y(t) is obtained as a Fourier transform 
(FT) of the signal.10 The power spectrum provides information on the amplitude of the 
predominant frequencies present in the time series. This information allows one to identify 
periodic, multi-periodic, quasiperiodic or nonperiodic signals. Usually the logarithmic power 
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spectrum plot is used to better distinguish between the broadband and periodic components. 
If the periods present are not of primary interest, then it is also used to evaluate the overall 
behavior of the time series. A power law dependence, and thus a linear dependence on a 
log- log plot, of S(f) given by 

s(f) ~ r p (i) 

that follows from the squared amplitude of the Fourier transform of the signal 



S(f) = lim 



/ - 



T 



e 2mft y(t)dt 



(2) 



is a hallmark of a self-affine phenomenon that is underlying the data. The properties of 
the signal can be further classified as persistent or anti-persistent fluctuations, depending 
on the values of the spectral exponent (3 in Eq. (1). If the signal possesses a tendency for 
repeating the sign of its fluctuations and, therefore, being of persistent type, then 2 < (3 < 3; 
when fluctuations having opposite signs follow one another, i.e. when 1 < (3 < 2, the signal 
is said to be anti-persistent .lB A signal with a spectral exponent obeying 1 < j3 < 3 is a 
nonstationary signal with stationary increments. In this case, the time series created by the 
increments of a nonstationary signal has a spectral exponent (3 within the interval [-1,1], 
and so is a stationary signalJa The latter property facilitates application of techniques of 
analysis to the increment signal. However the FT method, which produces second-order 
statistics, is insufficient to describe in full the signal scaling properties, because higher order 
moments may not be negligible. Strictly speaking, power spectrum analysis is suitable only 
for stationary time series. For more information on spectral time series analysis see [0-1.1 



Note that spectral analysis provides information on the scaling properties of the signal at 
high-frequency/short-time-scales, in contrast to the detrended fluctuation analysis method, 
which is reviewed next. 



B. DFA method and time dependence of the correlations 

A method that relaxes the requirement of stationarity of the investigated signal is the 
detrended fluctuation analysis (DFA) method.@ The DFA method is a tool used for sort- 
ing out lona range correlations in a nonstationary self-affine time series with stationary 
increments. The method has been used previously in the meteorological fieldiiHiJ It 
provides a simple quantitative parameter - the scaling exponent a, which is a signature of 
the correlation properties of the signal. Let the signal y(t) be defined between the beginning 
to and end of the observations t M , i.e. in [t ,t M ]. The DFA technique consists in dividing 
a time series y(t) of length iV into an integer number N/fi of nonoverlapping boxes (called 



also windows), each containing // pointstMa (fi = 4,5, ). The local trend z{ri) in each 
box is defined to be the ordinate of a linear least-squares fit of the data points in that box. 
The detrended fluctuation function F 2 (fi) is then calculated using: 

1 //V \ 

F 2 M = - E [y(n)-z(n)} 2 for k = 0, 1, 2, 1 . (3) 

f X n=kfi+l \ A* / 
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Averaging F 2 (n) over the N/fi intervals gives the mean-square fluctuations that are 
assumed to follow a power law 

< F 2 (//) > 1 / 2 ~/i°. (4) 

The DFA exponent a so obtained represents the correlation properties of the signal: 
a = 1/2 indicates that the changes in the values of a time series are random and, therefore, 
uncorrelated with each other, as in a Brownian random walk sequence. If a < 1/2 then the 
signal is anti-persistent (anti-correlated), and a > 1/2 indicate positive persistency (corre- 
lation) in the signal. It has been shown by Heneghan and McDarbyS that the relationship 
(3 = 1 + 2a holds true for stochastic processes, i.e. for fractional Brownian walks. 

The a-exponent value that holds true for a certain time interval called the scaling range, 
is a characteristic of the correlations in the fluctuations of a signal y(t). It is of interest, 
however, to test whether the correlations maintain the same properties in shorter intervals 
within [to, ijw] or whether they change with time, as it should be anticipated for nonstationary 
time series data. In order to probe the existence of so-called locally correlated and decorrelated 
sequences,0 one can construct an "observation box" with a certain width, u, place the box 
at the beginning of the data, calculate a for the data in that box, move the box by Av 
toward the right along the signal sequence, calculate a in that box, and so on through the 
Nth point of the available data. Each a-value is assigned to the end point of the box because 
all points in the box are needed in order to obtain a. A time-dependent a may be expected 
for t ranging from v to N. This approach is suitable for cloud base height data because it 
can be expected to reveal changes in the correlation dynamics of the clouds at various times 
for a given time lae v. If such a time dependence is found, then a multifractal approach is 
to be suggested^ In fact, it is expected that the usual fractal dimensionB D measuring 
the roughness of a signal^ is directly related to a through!!! D = 2 — a, whence D is 
time-dependent as well. 



C. Multifractal aspects 

The scaling behavior of a signal y(t) can change in a nonlinear fashion with the statistical 
moments, i.e. can be characterized by different scaling exponents. Such a signal is multi- 
affine and can be described through multifractal measures. E§c3 One ap_proach which we will 
follow here consists in studying the intermittency^] and the roughness" of the signal. The 
intermittency is quantified adopting the singular measure analysis The first step that this 
technique requires is to define a basic measure e(l; I) as 

™ = < = o,i,...,*-i 0) 



The notion of intermittency has no canonical dennitionc3 and covers a variety of phenomena. 
For us, intermittency consists in the existence of large and rare fluctuations with some structure 
localized in space and time for which a peculiar temporal behavior is found between periodic and 
chaotic regimes. 
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where Ay(l; I) = y{t i+ i) — y(ti) is the small-scale gradient field and <> denotes an average 
over the iV data points 

i N-l 

<Ay(l;l)>=-J2\Ay(l;l)\. (6) 
JV 1=0 

Next we define a series of ever more coarse-grained and ever shorter fields e(r; I) where 
< I < N — r and r — 1, 2, 4, 8 Thus the average measure in the interval [I; I + r] is 

e(r;0 = ~'l>(l;O- (7) 
The scaling properties of the generating function are then obtained for through the equation 



r 



-K{q) 



<e(r-iy>~^-} , g>0. (8) 
Using ii"(g)-functions one can defined the generalized multifractal dimension!^ 

!>(,) = (9) 

Note that at q = 1 the l'Hospital rule is used to obtain D(l). 

Further, the multi-affine properties of a time-dependent signal y(t) can also be described 
by the so-called qth order structure functions^ 

(\y(t t+r ) - y(U)\ q ) ~ r« ff W, z = 1, 2, . . . , N - r (10) 

with t = A, r = t i+r — tj. Whence and K(q) describe the multifractal scaling properties 
of nonlinear dynamic processes. In a monofractal case, H(q) and and K(q) take constant 
values. 



D. Hierarchy of exponents ^(7) 

The multi- affinity of y(t) means that one should use different scaling exponents H(q) in 
order to rescale a signal in various scaling ranges. In other words, this also implies that local 
scaling exponent 7 exists^ in order to characterize the local singularity of the signal. The 
density of points iV 7 (A) that have the same local scaling exponent is usually assumed0 to 
scale over the time span A (for any r) as 

iV 7 (A) ~ A~ h{7) . (11) 



From Ref. [ |32|] the following relations are found: 



7h> = ^ (12) 

h{ 1 {q)) = l + q 1 {q)-qH{q) (13) 

The h("y(q)) function^ is as naturally adapted to describe multi-affine signals as the multi- 
fractal spectrum for multifractal objects." 

Several of the methods reviewed in this section have been recently used in the meteoro- 
logical field .EHHJHiiiHii 
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III. DATA AND DATA ANALYSIS 



Next we apply the above methods to study the nonstationarity of cloud base height 
(CBH) data sets and the correlations in their fluctuations on all measured time scales; data 
are obtained at the Southern Great Plains site of the Atmospheric Radiation Measure- 
ment Program of the Department of Energy using a Belfort Laser Ceilometer (BLC) Model 
7013C.i3 The ceilometer is a self-contained, ground-based, optical, active, remote sensing 
instrument with the ability to detect and process several cloud-related parameters, among 
them cloud base height, cloud extinction coefficient and cloud layer depth. The ceilometer 
system detects clouds by transmitting pulses of infrared light vertically into the atmosphere 
and analyzing the returned signals backscattered by the atmosphere. The receiver tele- 
scope detects scattered light from clouds and precipitation.^ The ceilometer actively collects 
backscattered photons for about 5 seconds within every 30-second measurement period. The 
BLC measures the base height of the lowest cloud detected between 15 and 7350 m directly 
above mean ground level. 

The cloud base height signal measured on Sept. 23-25, 1997 is plotted in Fig.l. Data 
consists of iV = 7251 data points, as for the purpose of this analysis we consider the record 
from 5:15 UTC on Sept. 23, through 17:40 UTC on Sept. 25, 1997. It is well representative 
of similar events often occurring on a shorter time interval. 

The probability density function (PDF) for the cloud base height signal (data in Fig.l) 
is shown in Fig. 2. The values of the abscissa, Ay = y(t + At) — y(t) are unnormalized and 
therefore, given in meters. The cases At = 30 s and At = 30 min are displayed in Fig.2. 
They correspond to short and long time lags. Both PDFs are strongly non-Gaussian. The 
double pyramid triangular shape with a width growing with increased time lag is similar to 
that found in other meteorological cases.00 

We first check the stationarity of the data by applying a power spectral analysis. The 
power spectral density of the cloud base height data set is shown in Fig. 3. The upper 
curve represents the spectral density that is smoothed by applying a moving average over 
equidistant intervals on the log-scale. This procedure is anticipated to be better suited to 
identifying scaling properties of the spectrum.i It should be noted that the procedure may 
lead to a slight change in the slope of the linear fit that may to alter the scaling exponent 
slightly, as seen here from (3 = 1.46 ± 0.08 to (3 — 1.33 ± 0.06 for frequencies lower than 
1/15 min -1 , i.e. ~ 10~ 3 Hz. Nevertheless, because 1 < (3 < 3, we can conclude that, the 
cloud base height data are nonstationary with stationary increments. Futhermore, because 
1 < /3 < 2, the signal is anti-persistent. Note that the non-smoothed spectrum is more 
suitable for detecting characteristic frequencies and periods in the raw signal. The error 
bars in this paper are calculated according to standard techniques.0 

Similar values for the spectral exponent, e.g. /3=1.28±0.1 and /3=1.49±0.08, are ob- 
tained in Ref. [ P^J4TJ for cloud base height data measurements during the Atmospheric 



Stratocumulus Transition Experiment (ASTEX). ASTEX was designed to clarify the tran- 
sition from stratocumulus to trade cumulus cloudsi in the marine boundary layer in the 
region of the Azores Islands.! 1 ! Studies of more cases are needed to clarify whether the slight 
difference in the scaling exponents between stratus over land, which is the subject of this 
paper, and stratocumulus clouds over ocean is physically and statistically significant or not. 
The DFA leads to an a exponent equal to 0.18 ± 0.002 for time lags less than 55 min, 
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followed by an exponent 0.13 ± 0.001, as indicated in the inset of Fig.4. Notice that the 
Heneghan and McDarby@ relationship (3 = 1 + 2a holds well for the shorter time lag region. 
Incidently we stress that two spectral regimes are hardly seen from the above spectral power 
density graph in Fig. 3. 

The time dependence of the correlations is next discussed. Results are plotted in Fig. 4 
as a function of the cloud life time for v = 7 h and Au = 30 min. The v value is chosen 
such that the finite-size effects are avoided, while the Av value is somewhat arbitrarily 
chosen for an adequate display. Smaller values of Au lead to rougher curves representing 
the time-dependence of a, nevertheless without significantly altering the overall results. 

The generalized fractal dimensions D(q) are shown in Fig. 5 and seen to decrease with q. 
The multi-affinity of the CBH signal is also observed from the nonlinearity of the functions 
qH(q), 7(g) and K(q) in Figs. 5 and 6. The value of the if (g)-function at q = 1 defines the 
nonstationarity parameter Hi that is a measure of the roughness of the signal, while C\ at 
q=l 



c x -- dK{q) 



dq 



(14) 



usually defines the degree of intermittency of the signal .Hil@@'il 

The value obtained for the H (q)-i unction at q — 1, H(q — 1) = H\ — 0.21 ± 0.02 is close 
to the a-exponent of DFA method, a = 0.18±0.002 of this study and is similar to the results 
for the cloud base height data measured during ASTEX, for whichE^ a = 0.24±0.002 for June 
18, 1992 and a = 0.21±0.005 for June 15, 1992 and for which§ H(q = 1) = H x = 0.23±0.04 
for June 14, 1992 and H(q = 1) = H x = 0.21±0.03 for June 15, 1992. If the cloud base height 
data series was a strictly defined monofractal signal, then one can expect that a = H\% 

Note that q = 1 is a special case that corresponds to a monofractal behavior. Using 
l'Hospital's rule in Eq. (|T4|) one can obtain the information dimension of the CBH signal, if 
it scales monofractal 



dK(q) 

D(l) = l- 



dq 



(15) 

q=l 



The dashed line in Fig. 5 defines the monofractal case D(l) = 1 — C\ = 0.90. Having in 
mind that C\ is related to the mean of e(r; I) for a one-dimensional field, we note that the 
singularities that contribute most to the < e(r; I) > occur on a set with fractal dimension 
-D(l). In order to compensate statistically for their sparse spatial distribution, in sharp 
contrast to a Gaussian process, these extreme events are far more rare but far more intense. 
There is no intermittency at all if D(q) = 1, i.e. C\ = 0. 

The multi-affine properties of the signal are represented by the ^(7) functions following 
Eq. (|13D in Fig.7. Although the ^(7) -curve does not reach unity, it tends to some maximum 
value that corresponds to 70 ~ 0.21, a value equal to the a exponent of the signal (Fig. 4). 
The crossing of the x-axis by the /i(7)-curve defines the minimum value 7 m j n ~ that is 
related to the minimum value of the a exponent that is contained in the signalS Here the 
value of h(0) = 0.70, in agreement with the value l-qH(q) at its maximum as obtained from 
Fig.6, (see Eqs. (|T2p and (0)). 
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IV. CONCLUSIONS 



It is often asked what the numerical values of such above exponents mean, and whether 
they have any physical interpretation at all. Let us recall that the fractal dimension of a 
profile is a measure of its roughness, while generalized fractal dimensions have been given 
some thermodynamic interpretation elsewhereSS and here above (Sect. 2). The above a 
exponent indicates that the cloud bottom is quite rough. The occurrence of two scaling 
regimes might be understood either from a mathematical point of view or a physical one. In 
the former case, Hu et al.El have indicated that such an occurrence can be due to different 
noises and different trends. On the other hand, the presence of two scaling regimes well 
separated by a crossover time lag may indicate that two processes occur in stabilizing, or 
not, the cloud base height data set (or profile). Indeed we had observed that the a-exponent 
tends toward low values when the cloud is in quasi-equilibrium and reaches a high value 
when the cloud breaks apart0. We may conjecture here that the two scaling exponents 
describe two regimes of long time span (slow) correlations, which maintain stable droplets, 
and one at shorter time scales for fast correlations between droplets for both agglomeration 
or cloud fracture process. More cases need to be studied, however, in order to provide a 
more precise interpretation of these so close values of the a-exponents. 

Let us also recall that the simplest interpretation of a stochastic signal is through the no- 
tion of (fractional or not) Brownian motion, or random walk. In so doing, one can interpret 
the CBH signal as mimicking the distance from the origin (the initial time) traveled by a 
particle diffusing, e.g. on a lattice. An anomalous diffusion occurs if the signal correlations 
are not strictly characterized by the laws governing classical Brownian motion, but by other 
types of power laws. If those power laws are characterized by time-dependent exponents, 
then multifractality is expected. Usually the fractal dimension of a strange attractor mea- 
sures the minimum number of components of the phase space necessary to describe the 
dynamic process. No such fractal dimension has been studied here. It would be interesting 
to look for it in order to define how many interrelated physical parameters are necessary to 
describe the CBH variability, i.e. pressure, temperature, humidity, wind velocity, etc. Fol- 
lowing that line of thought, we can consider the CBH as representing stochastically driven 
turbulent eddies, in which phase transitions occur in three-dimensional space. In so doing 
the CBH is exactly a visual observation of these transitions in eddies moving up and down 
and horizontally in a stochastic way, as a random walk. In our case, the data set being 
studied is not really representing a particle in a 3D space but rather a projection on a ID 
space; only the bottom of the cloud, projected on a vertical axis is studied with respect to 
its up and down motion, - just like particle undergoing (fractional or not) Brownian motion 
in a plane. 

In summary, after recalling a few statistical analysis techniques, we have applied them 
to a typical continental stratus cloud base height profile data set. We have demonstrated 
that the CBH profile is a nonstationary anti-persistent signal with stationary increments. 
The spectral exponent found has a value similar to the one for stratoculumus clouds over 
the ocean reported by some of us in another study. The multi-affine scaling properties of 
the data series found reflect the complexity of the processes that produce them. Further 
work should be directed toward relating the scaling properties expressed through these 
statistical parameters to the dynamical properties of the clouds, an important step toward 
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understanding, modeling and predicting their dynamical behavior. 
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Figure Captions 



Figure 1 - Cloud base height signal measured on Sept. 23-25, 1997 at the Southern Great 
Plains site of the Atmospheric Radiation Measurement Program of the Department of En- 
ergy. The abscissa marks the time in hours after UTC on Sept. 23, 1997. The data series 
contains 7251 data points. 

Figure 2 - Probability density function P(Ay, At) (PDF) of the cloud base height signal 
measured on Sept. 23-25, 1997 (data in Fig. 1). The values on the abscissa, Ay, are given 
in meters. Triangles correspond to short time lag At = 30 s equal to the discretization step 
of the data series, and circles mark the long time lag PDF with At = 30 min. None of the 
PDF curves is displaced. 

Figure 3 - Power spectrum S(f) of the cloud base height signal measured on Sept. 23-25, 
1997. The upper curve represents the smoothed spectra, vertically displaced by two decades, 
that scales with a spectral exponent (3 = 1.33 ± 0.06. 

Figure 4 - The local a-exponent (black circles) given by the Detrended Fluctuation Analysis 
method as a function of time for the cloud base height (CBH) data measured on Sept. 23-25, 
1997. Error bars are indicated. The rescaled (divided by a factor 3000) CBH signal is shown 
by dots. Inset: Log-log plot of the DFA-function Eq.(|J), showing the scaling exponents 
a% = 0.18, c*2 = 0.13, and the crossover time lag at 55 min. 

Figure 5 - Hierarchy of generalized dimensions D(q) for the CBH data in Fig.l. The 
straight line is drawn to enhance the value q s at which the D(q) function starts to deviate 
from a linear dependence. The dashed line defines the monofractal case D(l) = 0.90; from 



Eq.(14) one thus finds C*i=0.10. The K(q)-f unction defined through Eq. (§) indicates the 



intermittency of the signal. 

Figure 6 - The hierarchy of exponents qH(q) and 7(g) indicating the multi-affine properties 
of the continental stratus cloud base height (CBH) data measured on Sept. 23-25, 1997. For 
q = 1, H\ = 0.21 ± 0.02 is a parameter of nonstationarity. 

Figure 7 - The /i(7)-curve for the stratus continental cloud base height data measured on 
Sept. 23-25, 1997; data in Fig. 1. 
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FIG. 1. Cloud base height signal measured on Sept. 23-25, 1997 at the Southern Great Plains 
site of the Atmospheric Radiation Measurement Program of the Department of Energy. The 
abscissa marks the time in hours after UTC on Sept. 23, 1997. The data series contains 7251 
data points. 
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FIG. 2. Probability density function P(Ay, At) (PDF) of the cloud base height signal measured 
on Sept. 23-25, 1997 (data in Fig. 1). The values on the abscissa, Ay, are given in meters. Triangles 
correspond to short time lag At = 30 s equal to the discretization step of the data series, and circles 
mark the long time lag PDF with At = 30 min. None of the PDF curves is displaced. 
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FIG. 3. Power spectrum S(f) of the cloud base height signal measured on Sept. 23-25, 1997. 
The upper curve represents the smoothed spectra, vertically displaced by two decades, that scales 
with a spectral exponent j3 = 1.33 ± 0.06. 
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FIG. 4. The local a-exponent (black circles) given by the Detrended Fluctuation Analysis 
method as a function of time for the cloud base height (CBH) data measured on Sept. 23-25, 
1997. Error bars are indicated. The rescaled (divided by a factor 3000) CBH signal is shown by 
dots. Inset: Log-log plot of the DFA-function Eq.(||), showing the scaling exponents a\ = 0.18, 
a 2 = 0.13, and the crossover time lag at 55 min. 
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FIG. 5. Hierarchy of generalized dimensions D(q) for the CBH data in Fig.l. The straight 
line is drawn to enhance the value q s at which the D(q) function starts to deviate from a linear 
dependence. The dashed line defines the monofractal case D(l) = 0.90; from Eq.(|l4|) one thus finds 
Ci=0.10. The /C(g i )-function defined through Eq. (@) indicates the intermittency of the signal. 
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FIG. 6. The hierarchy of exponents qH{q) and 7(g) indicating the multi-affine properties of 
the continental stratus cloud base height (CBH) data measured on Sept. 23-25, 1997. For q = 1, 
H\ = 0.21 ± 0.02 is a parameter of nonstationarity. 
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FIG. 7. The /i(7)-curve for cloud base height data measured on Sept. 23-25, 1997. 
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